Published in Journal of Physics G: Nuclear and Particle Physics. Copyright 1999 lOP Publishing Ltd 



Electromagnetic showers in a strong magnetic field 

V Anguelov and H Vankov 

Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 
72 Tzarigradsko Chaussee, 1784 Sofia, Bulgaria 

Abstract. We present the results concerning the main shower characteristics in a 
strong magnetic field obtained through shower simulation. The processes of magnetic 
bremsstrahlung and pair production were taken into account for values of the parameter 
X 3> 1. We compare our simulation results with a recently developed cascade theory 
in a strong magnetic field. 



1. Introduction 

Electromagnetic showers are a universal phenomenon. Besides occurring in matter or 
radiation field cascade, multiplication of electrons and photons can arise in a strong 
magnetic field. Such super-strong fields (~ 10^^ G) probably exist in the vicinity of some 
astrophysical objects such as pulsars, for example. In this case rotating neutron stars 
induce strong electric fields above the polar cap. Accelerated by these fields, high-energy 
particles (with energy up to ~ 10 TeV) move along curved magnetic field lines and emit 
curvature photons. The energy of these photons is enough to produce electron-positron 
pairs in magnetic and electric fields. The subsequent quantized synchrotron radiation 
by pairs will convert to a second generation of pairs and then an electromagnetic 
cascade develops in the pulsar magnetosphere. The shower development determines, to 
a considerable extent, the properties of the observed radiation from these objects. This 
was proposed for the first time in [|^ . Since electromagnetic cascades in strong magnetic 
fields were considered in many works mainly in connection with specific models of radio 
pulsars, gamma-ray bursts, blazars (see, for example, 

It is well known that the essentially non-zero probabilities for magnetic 
bremsstrahlung and pair production require both strong field and high energies 0]. 
The relevant parameter determining the criteria for this is: 

_ e H 
mc^ Her 

where e is the particle energy, H is the magnetic field strength, m is the electron mass 
and Her = 4.41 x lO^^ q 

The total probabilities (cross sections) for radiation and pair production for a given 
value of the magnetic field strength depend only on x and are shown in figure 1. One can 
see that magnetic pair production has significant probability for x > 0.1 (photon energy 
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Figure 1. The total cross sections for magnetic bremsstrahlung and pair production 
as function of x- 



must be > 2mc^). For effective shower development one needs even higher values of x 
(x ^ 1) because with increasing x the radiated photon spectrum becomes harder. For 
X > 1 (quantum region) the energy of the radiated photon is of the order of the electron 
energy. It is interesting to note that for a photon with energy ^ 7, 5 x 10^^ eV even 
Earth's magnetic field (~ 0.3 G) is strong enough {x ~ 1) to be a good environment for 
creating an electromagnetic shower. If such extremely high-energy photons are presented 
in the primary cosmic ray flux they will undergo cascading in the geomagnetic field before 
entering the Earth's atmosphere. This problem has been intensively discussed recently 
in connection with the detection of the highest cosmic ray events 0. 

As mentioned above, most of the treatments of cascades in a magnetic field were 
connected with specific models of astrophysical objects and are numerical in nature. 
The most general treatment of cascade properties emphasizing an analytical approach 
is made in |Q, where the steady-state kinetic equations for the electron-positron and 
photon distributions are solved in a strong magnetic field. 

A different approach to the shower study in a magnetic field is applied in 
It is similar to those for showers in matter and was motivated by the study of the 
primary gamma rays with extremely high energies (> 10^° eV) propagating through the 
geomagnetic field and Earth's atmosphere. The average shower characteristics obtained 
by numerically solving the system of cascade equations, show some of the main features 
of the cascade. While the shower is similar to those in matter for x ^ 1, its nature 
changes sharply for x < 1 which is connected with sharp increase of the photon free 
path. 

Recently, a kinetic theory of electromagnetic showers in a strong magnetic field has 
been developed in a similar to the cascade theory in matter, in approximation A 
Electromagnetic shower theories have been developed since 1937 following the works 
of Bhabha and Heitler and Carlson and Oppenheimer |^. Landau and Rumer 
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developed a complete theory in approximation A. Their work contains the formalism 
which is widely used in later shower theories. 

In the shower theory the mathematical description of the cascade process is based 
on the Boltzmann kinetic equation for particle flux density. The system of integro- 
differential equations of the one-dimensional Landau-Rumer theory is universal because 
it describes the shower development in any substance. Only the expressions for 
pair creation and bremsstrahlung probabilities per unit length are different. Using 
asymptotic forms of both processes for very high energies in a strong magnetic field, 
analytic formulae similar to those of standard cascade theory in approximation A for 
one-dimensional shower characteristics are obtained in ^j. But, as in matter, the kinetic 
equations were solved within certain approximations. 

In this work we give the results from the Monte Carlo simulation of the longitudinal 
development of electromagnetic showers in a strong magnetic field for x ^ 1. We 
present shower profiles for different ratios of the primary and threshold energies, Eq/E, 
and energy spectra of shower particles at different depths. We analyse the behaviour 
of the shower maximum with respect to Eq/E. We compare our modelled results with 
theoretical ones in order to estimate theoretical approximations and the range of their 
validity. 



2. The probability functions 

The main elementary processes leading to particle multiplication in a magnetic field 
are magnetic bremsstrahlung and magnetic pair production. The corresponding 



probabilities per unit length are P, ITT 
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where e and uj are the electron and photon energy and u = ui = . Parameter 

ch(t)p]-^ (z/t) dt is a modified Bessel 
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X was defined above. Here h = c = 1. Ky {z 
function known as MacDonald's function. Parameter v can have any real or complex 
values, here = | and |. For simplicity we assume that the electron (positron) is moving 
perpendicular to the magnetic field H. As already mentioned, the probabilities of both 
processes are essentially different from zero under the condition that the parameter 
X 3> 1 , which means that the particle energy e » ec, where Ec = mc^^. It is 
mentioned in |11| that this condition corresponds to the conditions u/x << 1 and 
■"i/x << 1; i-G- one can use in (|1]) the asymptotic form of (z) for z « 1, 
ir,(z«l) = |M(i)- + .... Then expressions may be simplified: 
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Using u = uj/e (correspondingly u = s/uj) we can 
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The total probabilities per unit length for bremsstrahlung and pair production are 
given by 



Wr (e) = [ n {e, iu)dLU = -^ [ 

Jo £3 Jo 

= 5.642^, 

£3 

Wp [u) = r 7 {u, e) de = 1.467^. (5) 

Jo LUS 

Unlike the Bethe-Heitler probability, 7f{e, u) does not contain an infrared divergence 
and because of this Wr (e) is finite. 

As mentioned earlier, the probabilities of both processes were obtained using the 
asymptotic form of the function K^, (z) for z << 1. This explains the behaviour of the 
total cross sections as a function of the particle energy of power (— |) (at fixed H) and 
in figure 1 this is the region of x ^ 1- Expression (^ coincides with the expression of 
the photon attenuation coefficient given in Erber's review when the asymptotic form 
of auxihary function T(x) for x >> 1? T{x) ~ 0.60%^ , is used. Here x = 
is the photon energy. 
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3. Simulation 



To investigate the shower characteristics in a strong magnetic field we developed our own 
Monte Carlo code. The probabilities @) were used to sample energies of the secondary 
particles - photon in bremsstrahlung and electron in pair creation. We constructed 
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Figure 2. Shower profiles for electron induced showers. Numbers refer to the ratio 

tables with cumulative distributions from probabilities (0) through a small step by u. 
The mean interaction length is the inverse value of the total probability, and (^). 

It was shown in |^ that the average ranges of both electron and photon, with 

1 

energy E'q in a magnetic field if, are of the same order of magnitude, ~ —2-. The 

1 

quantity L = which includes the matter parameters (magnetic field strength H) 
could play the role of a radiation length. In this case Eq is the energy of the particle 
initiating the shower. It is important to note, that unlike matter, here L depends on 
the energy Eq of the primary particle. L is a relatively small quantity. For example, 
in magnetic field H = 10^ G for Eq = 10^ GeV, L = 14.86 cm, for Eq = 10^ GeV, 
L = 148.6 cm. 

We considered the problem as one dimensional, i.e. we assumed that the shower is 
only developing in the direction of the primary particle entering the magnetic field at 
t = 0. The distance is measured in units L. All shower particles were followed down to 
some threshold energy E. 

4. Results and discussion 

We performed simulation for various sets of primary and threshold energies {Eq and E) 
and different magnetic field strengths H. 

Our results confirmed the theoretical prediction that the above-mentioned quantity, 
L, plays the role of a radiation length. Similar to the standard shower theory under 
approximation A, the longitudinal cascade development is independent of an absorber 
(magnetic field) when distances are measured in radiation lengths and the average 
behaviour of a shower is expressed by a function of Eq/E. 

Shower profiles for electron- initiated showers and different Eq/E are shown in figure 
2. The corresponding curves for photons are very close to those of electrons and because 
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Figure 3. Shower profiles for electron (solid curve) and photon (dashed curve) induced 
showers. 

of this they are not shown. When the primary particle is a photon (figure 3), the shower 
maximum is shifted by ~ | r.l. (where r.l. denotes radiation length) deeper, which is 
close to the difference between the mean interaction lengths of primary electron and 
photon. 

The differential energy spectra of shower particles at different stages of the shower 
development are shown in figure 4. The depth of 2.25 r.l. is near the shower maximum. 
The spectra can be described by power-law of energy, ~ E~^, with 5 increasing 
with the depth and approaching 1 for t — oo. 

As can be seen, the typical distance over which the shower develops is a few radiation 
lengths. The depth of the shower maximum increases very slowly with the increase of 
Eq/E approaching a limit. In very strong fields, e.g. in pulsars, L becomes very small 
(~ 10~^ cm for Eq = 10^ GeV) which means that the shower has the longitudinal 
spread of the same order. This confirms the theoretical prediction in [Q that the 
strong magnetic fields are effective screens for very high energy electrons, positrons 
and photons. 

It should be noted here that these properties of the shower are valid for such values 
of Eq, E and H that satisfy the condition x ^ 1 and the asymptotic expressions (2) 
can also be used. In addition, to consider the shower as one-dimensional, the electron 
energies must obey the condition ^ 

(^Y (—Y 

\mcV \H J 27r 

This criteria comes from the requirement that the electron gyroradius R = e/eH 
must be much greater than the typical length of the shower, i.e. L. 

However, detailed comparison of a numerical results shows that the Monte Carlo 
cascade curves differ significantly from the theoretical ones. This is illustrated in figure 
5 where shower profiles for the electron-induced showers and Eq/E = 1000 are shown. 
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Figure 4. Energy spectrum of electrons at three depths: line with squares-0.75 r.l., line 
with circles-2.25 r.l. (near the shower maximum), line with crosses-5 r.l. Eq/E = 1000. 
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Figure 5. Simulated shower profile (MC) compared with the theoretical curve in the 
adiabatic approximation (AA). 



This disagreement could probably be explained with the approximations used in 
to get a solution of cascade equations. In the case of a magnetic field, the shower 
theory is more complicated than the conventional theory in approximation A because 
the probabilities vr (e, tu) dtu and 'y{u!,e)de are not scaling functions, i.e. they do not 
depend only on the ratio of secondary to primary energies (functions of u/e or e/u 
in our notation). Mellin transforms lead to differential-difference equations for the 
distribution functions and its solutions are found by authors of [§] in the so-called 
adiabatic and modified adiabatic approximations. As pointed out in [^], the adiabatic 
approximation (AA) is limited to the region after the maximum. In the modified 
adiabatic approximation it is assumed that the explicit dependence of tt {e, uj) and 
7 (u;, e) on primary energy can be eliminated by its substitution with some mean energy 
per interaction in the shower. This leads to the same solutions with the modified 
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Figure 6. The number of electrons in tlie shower maximum as a function of 
Eq/E. MC-this work, AA-adiabatic approximation and BH-classical cascade theory 
in approximation A. 
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Figure 7. Depth of the maximum as a function of Eq/E. 

radiation length but, obviously, the shower development is distorted. If we use the 
probabilities 7r(£,a;) and 7(a;,£) in this modified way in our Monte Carlo code we 
obtain results very close to the theoretical ones. Figure 6 shows the number of electrons 
in the shower maximum A^max as function of Eq/E. One can see that the rise of A^max 
with Eq/E is slower than that of AA of the theory. A^max is exactly proportional to 
(Eq/E)^ . AA gives ^ close to one as it is in the standard theory with Bethe-Heitler 
cross sections (curve labelled BH in the figure) while our simulation (curve labelled MC) 
gives ~ §• Figure 7 shows the depth of the maximum tmax as a function of Eq/E for 
both theoretical and simulated showers. As can be seen, the behaviour of the modelled 
^max is too different from those of logarithmically increasing theoretical tmax- After some 
Eq/E the modelled tmax almost ceases to increase. 

It is easy to demonstrate this feature of showers in a magnetic field using Heitler's 
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Figure 8. Fluctuations of the number of electrons (a j (N)) as a function of the depth 
for different E/Eq. 
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Figure 9. Fluctuations of the number of electrons in the shower maximum {a/ {N)) 
as a function of E/Eq. Dotted line-showers in a magnetic field, solid line-BH showers. 



elementary cascade model |]I2|. In this model, electromagnetic particles subdivide into 
two particles with half the initial energy. In matter this takes place on every radiation 
length. In a magnetic field, however, the situation is substantially different. After each 
interaction the interaction length, which is a function of energy power of |, decreases due 
to particle energy splitting. If the primary particle has an energy Eq and the interaction 
length to then after the first interaction the number of particles will be two, each 
of energy Eq/2. The next interaction length will be and the number of particles 

four, each of energy -E'o/4. After interactions the particle number is already 2^ and 
their energy E = Eq.2~^ . The next interaction length will be — ^-l- The distance at 

H 

which this takes place is 
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The expression in the square brackets is a sum of a geometrical progression with 
q = 2~^^^ and thus 
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This expression shows that if N increases, i.e. Eq/E increases, tmax approaches a hmit. 
This means that after some N (or Eq/E) tmax practically ceases to increase. In our 
rather simplified model the maximum tma.x is 4.847to r.l. The first interaction length to 
is 0.117 r.l. for the electron and 0.682 r.l. for the photon. If we take the greater value 
this will lead to an estimation of the maximum t^ax of 3.31 r.l. whose value does not 
contradict with the curve modelled in figure 7. 

Another important characteristic of electromagnetic showers in a strong magnetic 
field which can be easily obtained from the simulation are fluctuations in the shower 
development. There is no theoretical consideration of this problem. The problem is 
complicated even for the standard cascade theory (see, e.g. fl^). 

In figure 8 the fluctuations of the number of shower electrons (cr/ (A^)) as a function 
of the depth for different Eq/E are shown. The primary particle is an electron. The 
behaviour of the curves is very similar to those of BH showers but the fluctuations in a 
magnetic field are significantly larger. This is illustrated in figure 9 where fluctuations 
in the shower maximum as a function of Eq/E are shown, compared with BH showers. 
Calculations for BH showers are performed by direct MC simulation in air for very 
high energies Eq and E, i.e. at the conditions where the standard cascade theory in 
approximation A is valid. 

The larger fluctuations in a magnetic field compared to BH showers are not an 
unexpected result. The main sources for the fluctuations of the number of particles in 
the cascade are fluctuations of interaction lengths and the random energy distribution of 
created secondary particles. Unlike the BH cross sections the mean interaction lengths 
for magnetic bremsstrahlung and pair-production processes are functions of the particle 
energy of power | (see (@J^)). The differential cross sections (^ depend on the particle 
energy as well. Their features lead to a significant probability that the secondary electron 
or photon will take a large fraction of the primary energy. 
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5. Conclusions 

A direct MC simulation of the longitudinal development of electromagnetic showers in a 
strong magnetic field has been performed. The processes of the magnetic bremsstrahlung 
and of the magnetic pair production were included in the simulation with asymptotic 
expressions for both probabilities valid for very high energies. Simulated results were 
compared with the theory of showers in a strong magnetic field developed in The 
main predictions of this theory - the dependence of the radiation length on magnetic 
field strength and the energy of the primary particle, the very small shower longitudinal 
spread and the closeness of electron and photon profiles, are confirmed by our simulation. 

However, the AA of the theory used in |^ was inadequate for a precise qualitative 
estimate of shower characteristics. 

Acknowledgment s 

The authors thanks A. Rekalo and T.Stanev for helpful discussions. This work was 
partially supported by a grant F-460 of the Bulgarian NFSR and by the Bulgarian 
Science and Culture Foundation. 

References 

[1] Sturrock P A 1971 Ap.J. 164 529 

[2] Daugherty J K and Harding A 1982 Ap.J. 252 337 

[3] Baring M 1989 Astron. Astrophys. 225 260 

[4] Bednarek W 1997 Mon. Not. R. Astron. Soc. 285 69 

[5] Erber T 1966 Rev. Mod. Phys. 38 626 

[6] Halzen F, Vazquez R, Stanev T and Vankov H P 1995 Astrophys. Phys. 3 151 
[7] Kanevsky B L and Goncharov A I 1989 Voprosy atomnoy nauki i techniki. Seria: Technika fiz. 
exp. 4 1 (in Russian) 

[8] Akhiezer A I, Merenkov N P and Rekalo A P 1994 J. Phys. G: Nucl. Part. Phys. 20 1499 
Akhiezer A I, Merenkov N P and Rekalo A P 1995 Nucl. Phys. 58 491 (in Russian) 

[9] Bhabha H J and Heitler W 1937 Proc. R. Soc. A 519 432 

Carlson J F and Oppenheimer J R 1937 Phys. Rev. 51 220 
[10] Landau L and Rumer G 1938 Proc. R. Soc. A 166 213 

[11] Bayer V H, Katkov B M and Fadin V S 1973 Radiation of Relativistic Electrons (Moscow: 

Atoniizdat) (in Russian) 
[12] Rossi B 1952 High Energy Particles (New York) 

[13] Uchaikin V V and Ryzhkov V V 1998 Stochastic Theory of High Energy Particle Transfer 
(Novosibirsk: Nauka) (in Russian) 



